APS/123-QED 



Open quantum systems approach to atomtronics 



O 
>v 



l> 



43 



> 

On 

\o 

O 

en 

in 

o 
o 



X 



R. A. Pcpino, J. Cooper, D. Mciscr, D. Z. Anderson, and M. J. Holland 

JILA, National Institute of Standards and Technology and Department of Physics, 

The University of Colorado, Boulder, Colorado 80309-0440, USA 

(Dated: May 19, 2010) 

We derive a quantum master equation to treat quantum systems interacting with multiple reser- 
voirs. The formalism is used to investigate atomic transport across a variety of lattice configurations. 
We demonstrate how the behavior of an electronic diode, a field-effect transistor, and a bipolar junc- 
tion transistor can be realized with neutral, ultracold atoms trapped in optical lattices. An analysis 
of the current fluctuations is provided for the case of the atomtronic diode. Finally, we show that it 
is possible to demonstrate AND logic gate behavior in an optical lattice. 
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I. INTRODUCTION 

The emerging field of atomtronics [1 3] aims to con- 
struct analogies of electronic components, systems and 
devices using ultracold atoms. In atomtronics, ultracold 
atoms move in an optical or magnetic potential in di- 
rect analogy with electrons moving in a semiconductor 
crystal. The motivation to construct and study atom- 
tronic analogs of electronic systems comes from several 
directions. 

First, the experimental atomtronic realizations 
promise to be extremely clean. Imperfections such as 
lattice- defects or phonons can be completely eliminated. 
This allows one to study an idealized system from which 
all inessential complications have been stripped. Conse- 
quently, one may obtain an improved understanding of 
the essential requirements that make certain electronic 
devices work. It is possible that a deeper understanding 
may feed back to the design of conventional electronic 
systems and could lead to future improvements. This 
lies parallel to the recent interest in single electron 
transistors in mesoscopic systems [4[ and molecules 
[5(, where many themes common with atomtronics 
emerge. A consequence of the near-ideal experimental 
conditions for optical lattice systems is that theoretical 
descriptions for atomtronic systems can be developed 
from first principles. This allows theorists to develop 
detailed models that can reliably predict the properties 
of devices. 

Second, atomtronics systems are richer than their elec- 
tronic counterparts because atoms possess more internal 
degrees of freedom than electrons. Atoms can be ei- 
ther bosons or fermions, and the interactions between 
atoms can be widely varied from short to long range and 
from strong to weak. This can lead to behavior that is 
qualitatively different to that of electronics jg-uj]. Con- 
sequently, one can study repulsive, attractive, or even 
non-interacting atoms in the same experimental setup. 
Additionally, current experimental techniques allow the 
detection of atoms with fast, state-resolved, near unit 
quantum efficiency [10(. Thus it is possible, in principle, 
to follow the dynamics of an atomtronic system in real 
time. 



Third, neutral atoms in optical lattices can be well- 
isolated from the environment, reducing decoherence. 
They combine a powerful means of state readout and 
preparation with methods for entangling atoms |ll|. Such 
systems have all the necessary ingredients to be the 
building blocks of quantum signal processors. The close 
analogies with electronic devices can serve as a guide in 
the search for new quantum information architectures, 
including novel types of quantum logic gates that are 
closely tied with the conventional architecture in elec- 
tronic computers. 



Fourth and finally, recent experiments studying trans- 
port properties of ultracold atoms in optical lattices |12| - 
Il4| can be discussed in the context of the atomtronics 
framework. In particular, one can model the short-time 
transport properties of an optical lattice with the open 
quantum system formalism discussed here. 



In this article we present a derivation of the master 
equation used to treat these specific open quantum sys- 
tems. Afterward, we provide a detailed analysis of atom- 
tronic analogies of the most elementary electronic com- 
ponents. These include conducting wires, diodes, and 
transistors. This work builds on a previous paper [3|, 
providing a comprehensive explanation of the underlying 
analytical and numerical methods, and additional anal- 
ysis of the components. Finally, we propose how AND 
logic gate behavior can be recovered in this open quan- 
tum system setting. 



This article is organized as follows. In Sec. [TT]we dis- 
cuss the general master equation formalism, introducing 
the specific systems to be investigated, the denning prop- 
erties of the reservoir, and the appropriate approxima- 
tions necessary to complete the derivation for the mod- 
els. In Sees. Mil and IIVI we apply the derived model 
to a variety of one-dimensional optical lattice systems in 
different open quantum system settings. The result is a 
collection of atomtronic devices that emulate electronic 
components. 
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FIG. 1: (Color online.) The illustration of an optical lat- 
tice, energetically tailored to have a transport response of an 
electronic diode, connected to two reservoirs whose respective 
chemical potentials induce an atomic transport from left to 
right. 



II. MODEL 

A. System 

Unlike the typical behavior of their electronic coun- 
terparts, atomtronic devices operate at the few-particle 
level. The necessary repulsive correlations between par- 
ticles are generally caused by either strong interactions 
for bosons, or quantum statistical effects due to the Pauli 
exclusion principle for fermions. In this article we focus 
on bosonic current carriers, partly to draw contrast to 
the electronic case. 

The optical lattice provides a clean, controllable envi- 
ronment for atomtronics components. Strong atom-atom 
interactions can be precisely tuned if the atoms are con- 
fined in a tight optical lattice. Holding the atoms in a 
lattice has two primary advantages. First, the strength 
of the interactions is enhanced due to the confinement 
within the lattice wells and second, by cooling the atoms 
deep into the lowest Bloch band regime, the center of 
mass motion of the atoms can be reduced to hopping be- 
tween neighboring lattice sites. This results in a simple 
theoretical description as well as a clean experimental 
realization. 

The dynamics of these ultracold atoms in this regime 
are accurately described by the Bose-Hubbard model 
[lol Uq |. In the lowest Bloch band, the Bose-Hubbard 
Hamiltonian is 

N 



Hbh = J2 ( ^ + \ U ^i (& - 1 



(id) V J 



(1) 



where i and j are lattice site indices, N is the total num 
ber of lattice sites, e, is the site energy, [/, is the on-site in- 
teraction energy, J^- is the hopping energy, the sum (i,f 
is taken between adjacent lattice sites, and N, 



Here, aj and a,i are bosonic creation and annihilation op- 
erators respectively for an atom in the lowest Wannier 
orbital centered at site i. Due to the precise tunability of 
the experimental parameters, e,, Ui, and J,j in the lab- 
oratory today, we use the Bose-Hubbard Hamiltonian to 
model our custom lattice atomtronic systems. 



B. Reservoir 

The reservoirs serve two purposes in atomtronics: they 
are 'sources' and 'sinks' of particles. The reservoirs them- 
selves could be experimentally realized in a variety of dif- 
ferent ways. For instance, a reservoir 'source' of atoms 
could be a f-dimcnsional optical lattice (Fig. [T]), or a 
strongly-interacting harmonically trapped ultracold gas. 
A reservoir 'sink' of atoms could be a coupling of the 
system to vacuum modes of atoms in an untrapped state 
or the densely spaced modes of a nearly empty potential 
well. For modeling purposes, suitable reservoirs must 
meet several requirements. The system-reservoir corre- 
lations must decorrelate from the system faster than the 
time scale over which the state of the system changes ap- 
preciably. This allows one to make the Markov approxi- 
mation. We will also assume that system and reservoirs 
are weakly coupled so that we can make the second-order 
Born approximation. In order for the reservoirs to serve 
as incoherent sources of atoms, their density of states 
must vary slowly over the spectrum of the system. 

In the model used in this work, it is assumed that 
the reservoirs are strongly-interacting bosons and at low 
temperatures, the states in the reservoir are filled to the 
chemical potential, and all states above the chemical po- 
tential are empty. This is analogous to the situation in 
a semiconductor crystal, where to a good approximation 
the electrons occupy all states up to the Fermi energy. 
Labeling the single particle excitations of the reservoir 
by a quantum number £;, we can write the Hamiltonian 
of the reservoir as 



Hr 



k 



T .kR k Rk, 
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where e k and R k are the energy and creation operator 
for the fc'th reservoir mode, respectively. 

Such a reservoir can be constructed out of bosons 
trapped in a deep optical lattice potential. If the in- 
teractions between the atoms are very strong the system 
enters the Mott-insulator regime and by adjusting the 
chemical potential one can achieve a situation were each 
lattice site is occupied by precisely one atom. The atoms 
in this system could then be coupled to the atomtronic 
device to furnish an incoherent source of particles, as de- 
picted in Fig. [TJ By arranging the distribution of the 
energies e k in the reservoir it is possible to achieve a sit- 
uation reminiscent of the Fermi-sea reservoir. We may 
then represent the reservoir by a collection of uncoupled 
harmonic oscillators. 



C. Elimination of the reservoir: master equation 
formulation 



where L = Lg + L^ + Ly. Noting that 
equation can be written as 



1. this 



The quantum master equation approach is often used 
in quan tum optics for describing open quantum systems 
p7H20| . In essence, it allows one to calculate observables 
associated with the evolution of a closed system without 
having to account for the free evolution of the reservoir. 
In this section we provide a detailed derivation of the 
master equation formulation for the system and reser- 
voir described above. The derivation involves deriving 
an equation of motion for the reduced density operator 
for the system. Here we construct this equation of mo- 
tion in a Liouville representation, since (as we discuss be- 
low) this form makes clear a way to go beyond the Born 
approximation, a necessity in this zero-temperature, un- 
conventional setting. 

The coupling of the system to the reservoir is by means 
of the exchange of particles between the system and the 
reservoir. The Hamiltonian for this interaction is 



H, 



^9kqR 



]flq 



H.c. 



(3) 
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where a q is the annihilation operator for a particle in 
a system lattice site q and gu q is the coupling matrix 
element between reservoir mode k and site q. 

The Hamiltonian for the system-reservoir interaction 
is given by 



H — Hs + Hr + H 



v- 



(4) 



where H$ = Hbh for our model. For each part of the 
Hamiltonian we introduce an operation L x , defined by its 
action on an arbitrary operator O by 



iL x Q = — [0, H x ], where 



e {S, R, V}. (5) 



We denote the density matrix of the system and reservoir 
with p. From the full density matrix, the reduced density 
matrices as and a R of system and reservoir are defined 
by tracing over the reservoir and system Hilbert spaces, 
respectively, 



<rs,R = Tzr,s\P]- 
We define the projection operator IP by 



(6) 



(7) 



and its compliment by Q — 1 — P. Under these def- 
initions, P and Q satisfy the usual projection operator 
relationships P 2 9 = P0, Q 2 = Q0, and PQ0 = 0. Us- 
ing the projection operators, the reduced density matrix 
for the system can be written as 



a s = Tr R [¥p}. 



(8) 



To find the equation of motion for as, we start from 
the evolution of the full density matrix p: 



dp 
~dt 



-iLp, 



(9) 



dt 



P- 



(10) 



Acting with P and Q separates this equation into the 
coupled equations 



dFp 
~df 



Wp + PLQ(Qp) 



dt 



The solution of equation (|T2|) is 



-lie 
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»p(t - r)dr. 



(11) 
(12) 

(13) 



where we have assumed that the system and reservoir 
initially uncorrelated, i.e. Qp(0) = 0. Using this result 
to eliminate Q/3 in Eq. (fTTj) yields 



—f- = -iPLPp - PLQ / e~ lUljllT QLP p(t - r)dr. (14) 
dt Jq 

Tracing over the reservoir Hilbert space leads to the equa- 
tion of motion for the reduced density operator of the 
system, 

^1 = -iL s as(t) + I "&(T)a s (t-T)dT, (15) 
where we have introduced the memory kernel 



®{t) 



-TV, 
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i(L s +L H +Qln 



Tyfffi 



(16) 



The first term on the right hand side of Eq. (fTS"|) is the 
free system evolution while the second term describes 
the irreversible contribution due to the system-reservoir 
interaction. 

Equation (|15p is the exact master equation under the 
condition that the system and reservoirs are initially un- 
correlated. It achieves our goal, in principle: all the dy- 
namics due to the coupling to the reservoirs is encapsu- 
lated in 3%(t). Once the memory kernel is known we can 
calculate the evolution of the reduced density matrix of 
the system without having to take into account the reser- 
voirs explicitly. However, it turns out to be impossible 
to solve for the memory kernel exactly in even the sim- 
plest of circumstances, so approximations must be made 
to continue from this point. 

Calculation of the memory kernel and its action on 
the system density matrix is straight forward under the 
Markov and Born approximations. The Markov approxi- 
mation assumes that the correlation time r c between sys- 
tem and reservoir is much shorter than the time scales 



over which the reduced density matrix of the system 
changes appreciably, i.e. 



da s ,. , 

Tc—TT < as- 

at 



(17) 



The Markov approximation consists of treating the sys- 
tem density matrix as a constant over time intervals of 
order r c , and accordingly we can pull it out of the integral 
in Eq. (fTS)) . The short correlation time also allows us to 
extend the limit of integration in Eq. (|15| to infinity. 

The Born approximation takes the memory kernel to 
second order in Ly and is thus a weak coupling approx- 
imation. In the conventional Born approximation, Ly 
is eliminated from the exponential term in the memory 
kernel. We cannot simply employ such an approximation 
here however since we are assuming that the strongly- 
interacting boson gas obeys a zero temperature Fermi- 
Dirac distribution characterized by fip. The hard edge 
at the Fermi energy causes logarithmic divergences in 
the second-order energy shifts of system levels as \xf ap- 
proaches system resonances. This divergence is due to 
the fact that, under the second-order approximation, as 
Hf approaches a system resonance, the number of reser- 
voir modes that the system is coupled to goes to one — not 
a continuum of modes. This would induce Rabi flopping 
of the atom between the system and the reservoir and 
not the irreversible reservoir action intended. In reality, 
the interaction of the system with the reservoir taken to 
all orders mixes the modes and leads to decay, in this 
situation. Exact simulations on small systems show that 
we can recover the proper dynamics at the hard edge 
by including the influence of higher order terms of the 
memory kernel expansion in Ly. We find that the fourth 
order term in the expansion provides a good estimate for 
the higher order corrections of the full memory kernel, 
yielding 

SHj) « -Ty r \L v L v (-T)a R e- [i< ^ s+hR)+r > ]T ] , (18) 

where Ly(r) = exp[— i(hg + L^)r]Lv r . The evaluation of 
the fourth order term shows that the decay r\ has a small 
dependence on the value \i R about the system eigencn- 
ergy difference. Here, we take r\ to be the mean value 
of the decay rate, which is where \x R equals the system 
eigenenergy difference, that is r\ = £3J 2 Tr/(2h 2 ), where 
@ is the density of reservoir modes, assumed to be con- 
stant in the region of interest. Introducing the rate at 
which particles enter or leave the system r = £>J 2 /h 2 , 



we have r\ = -nY^JI. Taking this value for r\ is reason- 
able because r\ is only important in the calculation for 
only small deviations. The modified memory kernel in 
Eq. (fT8"|) captures the correct long-time behavior of the 
exact memory kernel. 

With these approximations the master equation be- 
comes 



das 
dt 



/-•OC 

-iL s a s - / c~ r ' T Tr R [L v L v (-T)a R ]dTa s (t). 
Jo 

(19) 

Inserting the Liouvillians into Eq. (fT9| yields 



das 
dt 



ih s a s 
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dre-" T x 



Tr, 



Tr, 



d\R k , [a q {-T)Rl(-T),&sa R ] 



a q Rl, [al{-T)R k (-T),a s a R ] 



where we have used that 



Tr R [a R RlR k ,}=d kt k'(RlRk), 
Tr R [a R R k Rl,]=5 k , k ,(R k Rt), 



(20) 



(21) 



due to the fact that the reservoir is in thermal equilib- 
rium. 

We project Eq. ([20]) onto the energy eigenbasis of the 
system and trace out the reservoir degrees of freedom. 
This allows us to evaluate the r integral and to find a 
more explicit form of the master equation. Given two ar- 
bitrary system energy eigenstates \a) and \b) and adopt- 
ing the notation (a|G>|6) = (0) a &, we have 



{a\L s as\b} =uj a b(as)ab, 
(a\a' i (-T)\b) = exp(-iw ah r)(a t ) ab 



(22) 



with foujab = (E a — Eh) being the difference between the 
cigcncncrgies of | a) and 1 6) . Similarly, 



(z|E t (-r)|j)=cxp(-z^r)( J R t ) i 



(23) 



for the reservoir. Performing the integral over r produces 
the following closed form of the master equation: 



d(vs)ab 

dt 



-iuj ab (d- s ) ab - -^2 \ (ag)oc (ri. )cd(d q )cd {as)db— (rL )ac{a q )ac {os)cd{a\)db 

c.d ^ 

+ (vs)ac (r_)cd(a q )cd\(al)db~(al)ac(vs)cd (r )db(a q )db + (a q ) ac {?+) edited (ds)db 

- (F+ )ac(a q )ac (os)cd(a q )db + (vs)ac (1 + )cd(a q } c d (a q )db - (a q )ac(os)cd {?+ )db(a q )db H24) 



where 
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Note that 
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(RlRk), (25) 



(R k Rt). (26) 
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±UJk — OJab 



-ir6(±0Jk-Uab), 



where V indicates that integrals are to be interpreted in 
the Cauchy principal value sense. The real parts of the 
r± matrices then give decay rates that agree with the 
Fermi golden rule result while the imaginary parts give 
rise to level shifts in analogy with the Lamb shift in the 
hydrogen spectrum. 



D. Reservoir model 

The detailed physics of the reservoirs influence the evo- 
lution of the system through the coupling matrix ele- 
ments gk q , the occupation probabilities of the reservoir 
modes (R' k Rk), and the density of states of the reser- 
voir. The reservoir model used below assumes fermion- 
ized strongly-interacting bosons. Thus, (R* k Rk) = 1 for 
all modes below the chemical potential of the reservoir 



and (R k Rk) = for all modes above. The coupling of the 
reservoir modes to the system states is a slowly varying 
function of the mode energy. We model it by a constant 
coupling up to some high energy cut-off uj c in order to 
quench the ultra violet divergence that would otherwise 
arise, taking the form, 



\gk q \ 2 = \g q \ 2 6{uj c -uj k ), 



(27) 



with 6(y) the Heaviside step-function. The high energy 
cut-off lu c is much larger than any relevant frequency of 
the system, and it does not affect the system's dynamics. 



III. ATOMTRONICS APPLICATIONS 

In order to analyze the response of specific optical lat- 
tice configurations connected to reservoirs with different 
chemical potentials, we consider the steady-state solution 
of Eq. (j24|) and then proceed to solve for the matrix ele- 
ments (d's)ab- Once (<7s)ab is known, expectation values 
of atomic currents (J q ) into (and out of) a system site 
q can be calculated from reservoir's influence of the sys- 
tem population rates. For the population of each state, 
d{as}aa/dt we sum up net rates out of the system state 
| a) and then subtract the net rates into the state. We 
then sum over all of the system states \a) to obtain the 
net current transport on system site q: 



W = EE \(al)ac[(r ( ° ut) ) c d(a q ) c d] (as)da + (&s)ac [(if ut) )cd(al) cd ] (a q ) 

a c,d '- 

-{a\) ac {as)cd (r„ )da{a q )da - (T+ )ac(dl)ac (crs)cd(a q ) 

= y^,},(J)ad(vs)da 



da 



da 



a d 
Tr s I JqGS 
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FIG. 2: (Color online.) Illustration of the analogy between (a) 
the electronic circuit of a wire (with some inherent resistance) 
connected to a battery and (b) its atomtronic counterpart. 



as the current operator for site q projected onto the 
system eigenbasis. Using this convention, the sign of 
Trs[Jq<7s\ reveals whether the current flows into (— ) or 
out of (+) the system. 
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FIG. 3: (Color online.) Current of the atomtronic wire, in the 
case of weak coupling to the reservoirs, as a function of the 
left chemical potential. The current monotonically increases 
in quantized jumps as the on-site interaction energy to put 
on an additional particle is overcome. The parameters used 
to model the atomtronic wire in this simulation are e/U — 3, 
J/U = 3 x 1CT 2 , and HT /U = 1(T 6 . 



A. Atomtronics analogy of a simple circuit 

Here wc analyze the atomtronics counterpart to a sim- 
ple circuit of a battery connected to a resistive wire. As 
seen in Fig. [5J the analogy of a wire is an energetically- 
flat optical lattice, with uniform tunneling rates and in- 
teraction energies ( ej = e, Uj = U, and Jy = J for all 
neighboring sites). For this system, that is very weakly- 
coupled to the reservoirs, we calculate the atomic cur- 
rent as a function of chemical potential difference. This 
numerical experiment is carried out by initially setting 
both left and right chemical potentials (pl and /j,r) to 
zero. We raise /x^ so that an atomic transport is induced 
across the system from left to right, we compute the cur- 
rent ( Jr) out of the right side of the system shown in 
Fig. |3l The current increases with the chemical poten- 
tial difference, but in quantized jumps that correspond 
directly to the left chemical potential overcoming the 
on-site interaction energy needed to introduce a greater 
number of atoms onto the left site. A closer examina- 
tion of the numerical simulation implemented in Fig. [3] 
reveals two subtle features. Moving from left to right 
across graph in Fig. [3j the first current jump, occurring 
at (hl — e)/U — 0, the current increases in two steps 
where one might expect to observe a single jump, since 
the condition to put a particle in the lattice is /i£ > e. 
This is a result of the fact that the degeneracy of the 
Fock states in the one-particle manifold are split by 2J 
in the system state eigenbasis. In addition the jump in 
current is broadened slightly by the nonzero r\ and the 
system-reservoir coupling. This broadening, which makes 
the jump in current a smooth transition is more appar- 



ent in the atomtronic devices presented below where the 
system reservoir-coupling is taken to be orders of mag- 
nitude larger. Although the exact details for the second 
and third jumps are more complicated, the reasoning is 
the same: the eigenenergies are split by approximately 
2J, and the overall jump is smoothed out by r/. These 
are general properties of all of the numerical experiments 
described in this work. 



B. The atomtronic diode 

A diode is a device with an approximately unidirec- 
tional current characteristic. A voltage bias across the 
diode yields a current in one direction but not in the 
opposite direction if the voltage bias is reversed. Such 
behavior can be realized in an optical lattice by creating 
an energy shift in half of the lattice with respect to the 
other. We find that the diode characteristic persists as 
the number of lattice sites is increased. For simplicity, 
here we present the diode in a two-site lattice system. 
For the simulations in the rest of this paper, we assume 
e/U = 3, J/U = 3 x lO" 2 , T H/U = 1(T 2 . 

In the Fock basis for a two site system, there exist 
three states in the two-particle manifold: |20), |11), and 
|02), where \nm) refers to n particles on the left site and 
m particles on the right. The external energies of the 
two sites (ei and e 2 ) can be chosen so that the eigen- 
states |20) and |11) approximately degenerate, leaving 
both states far detuned from |02). This configuration of 
the site energies is given by €2 = (\ + U. We refer to this 
as the "resonance condition" . 



(a) 



(b) 



3c + 2(7 
|21>> 



|12>\ 



r\/ 

\-U 2e 

■0) | 

\7 



2e + ZU 



121) 

/o\/ 



112} 



\ 

2c + 3C7 



2t + [/ 2e + U I 02 ) 2€ + C7 2e + [7 l 02 > 

20) " ' |11>, / 20) -► |11> » / 

\/ \o/ \ ' 

01) —£ |01) 



1 10) 



\/ 



1 10) 



\ f 



|00) 



|00> 



FIG. 4: (Color online.) (a) Energy schematic of the reverse 
bias dynamics of the two site atomtronic diode. The red ar- 
rows represent the system transitions from an initially empty 
system. The gray arrows represent all other possible transi- 
tions. Regardless of which state the system starts in, it evolves 
almost entirely to the |02) state, (b) Energy schematic of the 
forward-bias dynamics of the two-site atomtronic diode. Blue 
and green arrows illustrate current-bearing cycles, while gray 
arrows represent all other possible transitions. 



Figure 2] illustrates how the resonance condition gener- 
ates reverse-bias and forward-bias behavior in a two-site 
optical lattice. As seen in Fig. [Ha), if one holds the 
left reservoir chemical potential at /i^ = and raises the 
right reservoir chemical potential fj,R, the system will un- 
dergo a transition from |00) to |01). The states |01) and 
1 10) are separated in energy by U. As a result, most of 
the population remains in the 1 01} state. Increasing ijlr 
above the point where the transition from 1 01) to 1 02) is 
allowed, the system remains almost completely settled in 
the 1 02) state. 

As seen in the Fig. Qlb), if °nc holds \ir = and 
raises ^l, the system first undergoes a transition from 
1 00) to 1 10). However, increasing fi^ so the system evolves 
to 1 20) leads to a very different situation than in the 
above case: since |20) is resonant with |11), both states 
are simultaneously populated. Since /Ur = takes all 
particles out of the site on the right, the system can make 
a transition from |11) back to 1 10) . The combined effect 
of setting /xl and fiR to these values is to force the system 
to undergo a closed cycle of transitions between 1 10} , |20) 
and 1 11). The result is a net atomic transport (or current 
flow) across the system. A second contributor to the net 
current through the system is the fact that fiL allows 
transitions from |11) to |21). Thus, an additional current- 
generating cycle exists: |11) to 1 21) to |20) and back to 
1 11). Both cycles contribute positively to a net current 
flow across the system. 

For systems consisting of N lattice sites, the diode 
configuration consists of two connected, energetically-flat 
lattices whose energy separation is Ae = U. The dynam- 
ics do not change for larger systems since the degeneracy 
of the flat lattice allows for effective transport across the 
lattice, allowing a particle to enter the left site of the 
half of the lattice is energetically degenerate with find- 
ing a particle at the junction. Thus, there exist current 
cycles initially generated from the 2 2 2 ... 2 ) 0. ..()()()) to 
|222...1) (g) 1 1.. .000) transition. Going the other direction, 
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FIG. 5: (Color online.) 
lattice atomtronic diode. 



Schematic of the four-site optical 
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FIG. 6: (Color online.) Current responses of (a) forward-bias 
two site diode, (c) the four site diode as a function of jil 
(with jj,r — 0); current response of (b) the reverse- bias two 
site diode and (d) four site diode as a function of fj,R (with 
Mi =0). The chemical potentials of all figures are normalized 
to the resonance condition, and the currents normalized to To. 



one can go to |0...000) <g> |222...2), but conditions are not 
energetically favorable to allow atomic transport across 
the junction. Figure [5] is a schematic of a four-site atom- 
tronic diode in the forward-biased direction. 

Supporting the behavior that the dynamics of the 
diode are qualitatively-independent of the overall size of 
the lattice, Figs. [6^a-b) and EJc-d) are numerical simu- 
lations of the current responses of the two site, and four 
site diodes, respectively. The general features of both 
diodes are qualitatively identical. 

Given the current characteristics of the diode, one 
might ask how much does the signal fluctuate. Recall- 
ing the current operator from Eq. (|31|) . we can calculate 
the autocorrelation function for the current, 



(J q (r + t)J q {t)) = Tr s [J q ( T + t)j q {t)a s ], 



(32) 



using the quantum regression theorem. 

To simulate an actual measurement, we convolve this 
correlation function with an exponentially decaying filter, 
exp(— t/T). The Fourier transform of this convolution 
yields a time-averaged spectral density function S(u>,T). 
Our time-averaged signal-to-noise ratio SNR as a func- 











fP" 




e + U 




■// 




e + U+lftJ 


M 




e + U+J 

e + (/ + 3/2J 


l 




e + U + 2J 


1 






J 





(a) 



3, 



0.5- 

0- 

0.7 0.8 0.9 1 1.1 1.2 

FIG. 7: (Color online.) Current response of the atomtronic 
FET as a function of the chemical potential: small changes in 
£2 lead to appreciable changes in the current response. That 
is, as the system is detuned from its resonance condition, there 
is a fall off of the net current value across the device. 



tion of T is then 

SNR(T) 



(J q ) 



\j I(T ' S(uj,T)du 



(33) 



For long averaging times T we find, for the conditions of 

Fig. m 



SNR(T) 



(34) 



For typical optical lattice experiments, U > ■> 1kHz is 
achievable, which implies that \/8Tq ~ lOyHz. There- 
fore, a signal-to-noise ratio of 10 can be achieved by av- 
eraging the atomic current for about 1 second. 



C. The atomtronic field-effect transistor 

A field-effect transistor (FET) is a device that allows 
an externally-applied field to affect the current through 
the device. This characteristic allows the FET to be uti- 
lized as an amplifier. Since the diode is optimized when 
the resonance condition is imposed on the optical lat- 
tice, small deviations from the resonance condition lead 
to large changes in maximum current propagating across 
the lattice. This is precisely the behavior of a FET where 
a current is controlled by an applied voltage. In Fig. [3 
we plot several current results for the forward-bias con- 
figuration as the separation in the external energy of the 
second site is raised past the resonance condition by frac- 
tions of J, the smallest system parameter in the model. 



D. The atomtronic bipolar junction transistor 

A bipolar junction transistor (BJT) is a three terminal 
device in which the overall current across the emitter and 
collector is controlled by a much weaker current via the 
base. Two practical applications of the BJT arc signal 
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FIG. 8: (Color online.) (a) An illustration of the BJT lattice- 
reservoir system, (b) The energy schematic of the three-site 
optical lattice under the extended resonance condition. Here 
the left reservoir is set to maintain an occupancy of one atom 
on the left site and the right is set to remove all atoms. If 
the middle reservoir is set to remove all particles, then the 
system evolves to the |100) state (red arrow). If the middle 
reservoir is set to maintain an occupancy of one atom on the 
middle site, then the degeneracy between states 1 110} , |020), 
and 1 011) allows current to traverse the system. 



amplification and switching (on and off) of the emitter 
current. 

Realization of B JT-likc behavior in atomtronic systems 
requires at least three sites connected to three different 
reservoirs. If the atomtronic diode is considered an atom- 
tronic p-n junction, one might guess that the atomtronic 
n-p-n transistor would entail raising the external energy 
of the left and right (collector and emitter) sites higher 
than the middle (base) site by the on-site interaction en- 
ergy. This configuration is illustrated in Fig. [5Ja). The 
reason this configuration yields BJT-likc behavior is due 
to the approximate degeneracy between the Fock states 
|110), |020), and |011). 

We implement numerical simulations of this lattice 
configuration by fixing a chemical potential difference 
across the lattice. The left reservoir chemical potential 
Hl is set to maintain one particle on the left site and the 
chemical potential of the right reservoir /z/j is set to zero. 
The middle chemical potential hm starts at zero and is 
increased to allow a single atom to enter the middle site. 

When there are no atoms on the middle site, the con- 
figuration of the reservoirs pumps the system into the 
1 100) Fock state (as seen in Fig. Mb)). The |100) and 
the 1 001) states are degenerate with each other, but the 
system must undergo a second-order, off-resonant tran- 
sition via the |010) state from 1 100) to |001). Such tran- 
sitions are suppressed by a factor of (J/U) 2 and thus 
become less likely as the energy difference between 1 100) 
and 1 010) increases. Thus, when the middle reservoir is 
set to maintain zero atoms on the middle site, the net 
current out of the emitter is minimal. 

When the middle reservoir's chemical potential is in- 
creased to allow a single atom into the middle site of 
the system, the degeneracy between the 1 110) , |020), and 
|011) states is accessed, which allows atoms to travel 




FIG. 9: (Color online.) Characteristics of the atomtronic 
BJT. (a) For a fixed collector-emitter voltage bias, the cur- 
rent response measured through the base leg (green) and the 
emitter leg (blue) are plotted as a function of the base chem- 
ical potential, (b) The current out of the emitter is plotted 
vs. the current out of the base. Here a large linear gain is 
observed. 



across the system. One issue with this configuration is 
the following: in order to get to |011), the system has 
to make a transition through the |020) state. Since the 
middle reservoir is set to maintain an occupancy of one, 
but not two, atoms on the middle site, one of the atoms 
can be lost to the middle reservoir, leading to a loss of 
current out of the emitter. If the couplings of all three 
reservoirs to the system are equal, then the result is the 
current measured passing through the base turns out to 
be even greater than the current measured out of the 
emitter. Thus, the system represents an inefficient tran- 
sistor realization. On the other hand, if the middle reser- 
voir were to be coupled weakly compared to the other 
reservoirs, then current predominantly leaves the system 
via the emitter, which is the desired behavior. 

Figure |9ja) is a numerical simulation of the current 
out of the emitter and the base as a function of the base 
chemical potential. The coupling strength of the base 
connected to the reservoir is one fifth the collector and 
emitter reservoir coupling. It should be noted that the 
region where the proposed atomtronic transistor mimics 
the electronic BJT is limited to the transition region, 
or current jump. One can increase the length of this 
region by increasing the overall system-reservoir coupling 
strength. Figure [9jb) shows that the gain of this device 
is fairly linear. 



IV. DISCRETE ATOMTRONICS LOGIC 

Integrated circuits arc designed with a very large num- 
ber of transistor elements to perform a desired function. 
The demonstrated ability to realize atomtronic diodes 
and transistors thus motivates the question as to whether 
higher functionality can be realized with these ultracold 
atomic systems. Here we look at the most fundamental 
of these, the atomtronic AND logic gate. 

A traditional logic element is a device with a given 
number of inputs and outputs, composed of switches, 
that generates a series of logical responses. Such logi- 
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FIG. 10: (a) Electronic schematic for the AND logic gate, a 
device constructed by cascading two transistors in series, (b) 
Atomtronics AND schematic. The atomtronic AND logic gate 
is constructed exactly like its electronic counterpart, i.e. by 
cascading two transistors in series. A chemical potential bias 
I^l — fJ-R across the device attempts to drive a current across 
the device. No substantial current is observed, however, un- 
less both \xa and Hb supply atoms onto the base terminals of 
their respective transistors. 



cal behavior can be expressed in a truth table composed 
of l's and O's ('ons' and 'offs'). Logic elements are the 
fundamental building blocks of computing and discrete 
electronics. In table [IJa) the truth table for the AND 
logic gate is given as an example. The next level of com- 
plexity in emulating electronic systems is to create logic 
elements from the atomtronic components. 

An AND gate is a device with two inputs (A and B), 
and one output (O). As illustrated in table Uja), the 
device characteristic of the AND gate is that O remains 
off unless both A and B are on. In electronics, such a 
device can be constructed by connecting two transistors 
in series (as illustrated in Fig. fTUTa)). By analogy, if 
the atomtronic BJTs are connected in the same series 
configuration (as illustrated in Fig. flOfb)). the AND 
gate truth table can be generated. 

When constructing practical logic circuits, the values 
of the Is and 0s are not strict values, they are defined 
within a given range. The data in table Htb) has been 
generated in a numerical experiment of the configuration 
depicted in Fig. |T0](b). For this particular experiment, 
the maximal current out of the device is at least a factor 
of 6 greater than any other measured current out. Thus 
a discernible difference between 'on' and 'off' is observed 
and the output currents reproduce the AND gate truth 
table. Such a difference can also be enhanced by increas- 
ing the on-site interaction energy U of the lattice. 



V. CONCLUSIONS 

In this article, we have derived a general model for 
treating a specific class of open quantum systems where 
the reservoirs act as sources and sinks for particles mov- 
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Atomtronic AND Gate Simulation 
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0.00 


3.2 





0.01 





3.2 


0.16 


3.2 


3.2 


1.00 



TABLE I: AND gate truth table where the two inputs are 
labeled by A and B, the output labeled by O, and the 1 or 
entries represent 'on' and 'off', (a) Truth table describing the 
device characteristics, (b) Numerical simulation of the atom- 
tronic AND gate normalized to the maximal output current 
'O'. Here, the inputs A and B are the respective chemical po- 
tentials measured with respect to the onsite interaction energy 
(Ha/U hb/U) and the output is the atomic current normal- 
ized by the maximum output current (I/I ma x)- ln practice, 
the logical 'ins' and 'outs' are never really '0' or '1' but are 
defined by setting threshold values. Since the maximal output 
current is larger than all other output values by at least a fac- 
tor of 6, we can conclude that there exists logic gate behavior 
in this optical lattice setup. 



ing into and out of the system. Such a formalism can be 
used to study atomic transport across arbitrary multiple 
potential well configurations. Here, the formalism was 
used to show how neutral atoms in custom optical lat- 
tices can exhibit electronic diode, FET, BJT, and AND 
gate behavior. 

Looking forward, we aim to develop more compli- 
cated atomtronic devices such as additional logic ele- 
ments, flip-flops, and constant current sources by cascad- 
ing our current atomtronic components in a manner anal- 
ogous to the development more sophisticated electronic 
devices. The simulation of the AND gate is promising 
since it demonstrates the possibility of cascading atom- 
tronic components to make more sophisticated devices. 
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